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Abstract 

Contributions to B — ► ixlv decay from 1/mq order corrections are analyzed 
in the heavy quark effective field theory (HQEFT) of QCD. Transition wave 
functions of l/mq order are calculated through light cone sum rule (LCSR) 
within the HQEFT framework. The results are compared with those from 
other approaches. 

PACS numbers: ll.55.Hx, 12.39.Hg, 13.20.Fc, 13.20.He 

Keywords: B — * ttIu, 1/mg correction, effective field theory, light cone sum 
rule 



1 



Typeset using REVTfjX 



I. INTRODUCTION 



Heavy to light semileptonic decay B — > itlv has attracted much interest due to its spe- 
cial role in determining the important Cabibbo-Kabayashi-Maskawa (CKM) matrix element 
\V ub \. A lot of work concerning this decay may be found in the literatures [1-3]. Since B 
meson consists of a heavy quark and a light quark, it is natural to explore B decays with 
considering the heavy quark symmetry. This has been proved to be powerful in relating var- 
ious heavy hadron processes, and distinguishing between long and short distance dynamics 
reliably. 

The leading order wave functions of the decay B — > irlv have been calculated in Ref. 
[4] by using light cone sum rules within the heavy quark effective field theory (HQEFT) 
[5-9]. In this paper, we present the next to leading order study for B — > ttIv decay in 
HQEFT. In section II we formulate this decay up to 1/itiq order in HQEFT. In section 
III the I/ttiq order wave functions are calculated by using light cone sum rule method in 
HQEFT framework. Section IV devotes to the numerical analysis and relevant discussions 
in comparison with other approaches. A short summary is outlined in section V. 



II. B -► ttLu DECAY IN HQEFT 

The semileptonic decay B — > nlu is determined by the matrix element 
(n(p)\uYb\B(p B )} = 2/ + (gV + (U(q 2 ) + f-tfW 



= U(q 2 ) Pb + p"- " 2 V + fo(q 2 ) B 2 V (2.1) 
l qz \ qz 



with 



/o(g 2 ) = -T^/-(? 2 ) + / + (g 2 ), (2.2) 

where q = p B — p is the momentum carried by the lepton pair. 

In HQEFT, the heavy quark field Q is decomposed into two parts Q + and Q~ , which are 
formally corresponding to the two solutions of the Dirac equation. For the free particle case, 
they are known to be the quark and antiquark fields respectively. In order to shift out the 
large component of heavy quark momentum and also to make the leading order Lagrangian 
explicitly ijiq independent, an effective heavy quark field can be conveniently defined as 

Q+ = e iM Q v-XQ+ = e iM Q v-xp + Q+ (2.3) 

with v being an arbitrary four-velocity satisfying v 2 = 1, and P + = (1 + f)j1 being a 
projection operator. Under this definition, Q+ is actually the "large component" of heavy 
quark field. Similarly, one can define the "large component" of heavy antiquark field. This 
becomes transparent if one considers a free quark field in the frame v = (1,0,0,0), where 
one has in the momentum space 
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n+ 1 + i> ( \ E + m 1 1 ^ 



a 1 -i , x E + m ( \ . > 

Here w s and v s are the four- component spinors, while (p s and x« are the two-component Pauli 
spinor fields that annihilate a heavy quark and create a heavy antiquark, respectively. A 
detailed discussion on the "large" and "small" components of the heavy quark and antiquark 
fields can be found in Refs. [5,6]. 

Integrating out the small component of quark field as well as the antiquark field, the 
effective Lagrangian can be represented in terms of as 



c { ; f f = q + v iv ■ DQt + —q:(h? ± ) 2 q: + o{\). (2.5) 

itlq m,Q 



At the same time, a heavy-light quark current qYQ with Y being arbitrary Dirac matrices 
is expanded into 

qYQ - e—^{gYQt + -i-gT^^)^ + 0(-L)}, ( 2 .6) 



2m Q iJp\\ v m Q 

where q denotes an arbitrary light quark. In Eqs.(2.5) and (2.6) the operators ilp\\ and Hp±_ 
are defined by 

iIp\\=tfvD, ijp^ = iJp-iiP\ V (2.7) 

For detailed derivation of HQEFT, the interested readers are refered to Refs. [5,6]. 

For the B — > ir weak transition matrix element, it receives I/tuq order contributions not 
only from the expansion of effective current (2.6) but also from the insertion of the effective 
Lagrangian (2.5). As a result, the heavy quark expansion (HQE) of the transition matrix 
element can be simply written as [6] 



(n\uYb\B) = J^{(«\u^Q:\M v ) + -J—(v\u^-^P + (Dl 

V B Q 

+ l -a aP F^)Q + v \M v ) + 0(l/m 2 Q )}, (2.8) 

where A# is defined by the mass difference of B meson and bottom quark, A# = mg — m^, 
and F a/3 is the gluon field strength tensor. \M V ) is an effective heavy meson state defined in 
HQEFT to manifest the spin-flavor symmetry, 

{M v \Qt^Qt\M v ) = 2Av» (2.9) 

with the binding energy A = lim^^oo Am being heavy flavor independent. 

Note that in Eqs.(2.6) and (2.8) we have used the operator l/(iv ■ D) to effectively 
represent the contraction of effective heavy quark fields (or say heavy quark propagator) [6,9]. 
In calculating Green functions in the next section, we will treat this operator as contraction 
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of and Q+. It is seen in Eq.(2.8) that the 1/toq order corrections to B — > 7r transition 
are simply attributed to one kinematic operator and one chromomagnetic operator. 

Based on heavy quark spin-flavor symmetry, we parameterize the matrix elements in 
HQEFT as 



(n(p)\uTQ+\M v ) = -Tr[7r(v,p)TM v }, 



{n(p)\uT 



iv ■ D 



(n(p)\uT—^-a a(3 F^Qt\M v ) = -Tr[itf {v,p)TP + -<j aP M v l 



(2.10) 
(2.11) 

(2.12) 



where M. v is the flavor independent spin wave function for pseudoscalar heavy mesons, 

M v = -v / AP +7 5 . 

The functions ir(v,p), ni(v,p) and ^^{v.p) can be generally decomposed into 



ir(v, p) = 7 5 [A(v ■ p) + f>B(v ■ p)] , 

Tn(v,p) = l 5 [fa(v-p)+pfb(v-p)}, 



(2.13) 
(2.14) 

nf(y,p) = 7 5 [(pV - P P l a )(si(v ■ p) + f>s 2 (v ■ p)) + ta a ?(s 3 (v ■ p) + f>s A {v ■ p))\ (2.15) 

with p^ = p^/(v ■ p). The Lorentz scalar functions A, B, /„, f b and = 1,2,3,4) 
are independent of the heavy quark mass. Their scale dependence is suppressed in these 
formulae. In Eqs.(2.13)-(2.15) we use for convenience v -p as the variable of wave functions, 
which is related to the momentum transfer q 2 by 



f = v -p 



m 2 B + m 2 — q 2 



2m i 



Now after trace calculation Eq.(2.8) yields 



mWta,)) = 2^ '^(A + J-^) + + J- Bl ) + O(-L)} (2.16) 



with 



Ai = /„ + 2si - 2s; 



3s 3 , 



Bi = f b - 2si + 2s 2 - 3s 4 , 



(2.17) 
(2.18) 



where = m 2 is used. Eqs.(2.1) and (2.16) yield relations between the form factors f± and 
the universal wave functions, 



f±k 2 ) 



As a comparison, in the usual heavy quark effective theory (HQET) one uses the effective 
Lagrangian 
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C eff \ H QET = Q> ■ DQt + ^-Qt(il?±) 2 Qt + 0(±-), (2.20) 
and the heavy-light current expansion (at tree level) 

qTQ + - qTQ+ + J—qTi]p ± Q+ + 0(\). (2.21) 
2m Q v m Q ^ 

In some references the notation h v is used for the field variable here. 

For the heavy quark expansion of the transition matrix element in HQET, the I/tuq 
order corrections from the current expansion (2.21) and those from the insertion of the 
Lagrangian (2.20) can not be attributed to the same set of operators. In particular, there 
are both operators with even and odd powers of ilfi± which contribute at I/ttlq order. 
Therefore the corrections from these two sources have to be characterized by different sets 
of wave functions and considered separately. 

In HQEFT, however, there are only operators with even power of ilp± appearing in the 
effective Lagrangian (2.5) and current (2.6). The terms with odd powers of ilp± are canceled 
due to the inclusion of the antiquark contributions. Therefore the I/tuq order corrections 
to transition matrix elements from the two sources can be attributed to the same set of 
operators, and the final formulation in HQEFT is simple [6]. In this paper, we need only 
calculate the two composite functions A 1 and B ± . 



III. LIGHT CONE SUM RULES IN HQEFT 

To calculate the 1/mQ order wave functions Ai and Bi, we consider the two-point cor- 
relation function 

F" = iJ d 4 xe^- m ^ x (n(p)\T{u(x)r-^P + (Di + l -o a pF^)Q+{x), 

g+(0)z 7 5 rf(0)}|0) (3.1) 

where Q^(0)i7 5 <i(0) is the interpolating current for B meson. Inserting between the two 
currents in Eq.(3.1) a complete set of intermediate states with the B meson quantum number, 
one gets 

AxvP + Bip* i f°° j n p a (v ■ p,s)v^ + p b (v ■ p,s)p^ 



2Ar - uj 



+ r° ds Mvp,s^ + Pb[ v.p,s^ + subtraction (32) 

B — UJ J so S — UJ 

with u = 2v ■ k, where k = pb — tkqv is the residual momentum. F is the decay constant of 
B meson at the leading order of 1/rriQ expansion, defined by [8] 

(0\qTQ+\B v ) = jTr[TM v ]. (3.3) 

At the same time, in deep Euclidean region the correlator (3.1) can be calculated in effective 
field theory. The result can be written as an integral over a theoretic spectral density, 

s — UJ 
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In the sum rule analysis, one assumes the quark-hadron duality, and equals the hadronic 
representation (3.2) and the theoretic one (3.4), which provides a sum rule equation. Fur- 
thermore, in order to enhance the importance of the ground state contribution, to suppress 
higher order nonperturbative contributions and also to remove the subtraction, a Borel 
transformation 

b<->. lim t?ri ( * r 

-w.n^oo n\ dtu 

-u/n = T 

should be performed to both sides of the sum rule equation. With using the formulae 

B^_L- = e s/T^ B ( ^e^ = 5(X-^), (3.5) 

one gets 

2iF{A 1 tf + B 1 p»)e- 2KB ' T = r dse~ s / T [p a (£, s)u" + p b (£, s)p^} (3.6) 

Jo L J 

with 

Pa&sW + ptfos)? = B^B^F^u). (3.7) 

The operator l/(iv ■ D) in Eq.(3.1) may be effectively regarded as a heavy quark prop- 
agator [6,9]. Namely, in the following we conveniently calculate the correlator (3.1) via the 
three-point one: 

i J d 4 xe^- m Q v > x (n(p)\T{u(x)rQ:(x),i 3 J d A yQ+{y)(p\ + ^o- a pF a ^Q+(y), 

Qt(0)tl 5 d(0)}\0). (3.8) 

In calculating the three-point function (3.8), we represent the nonperturbative contributions 
embeded in the hadronic matrix element in terms of light cone wave functions up to twist 
4. Among them are the two-particle distribution functions defined by [1,2,10] 

(7^)1^)7^(0)10) = -ijTU C due iu ^[^{u) + x 2 9l (u)] 

Jo 

) / due iu ^g 2 (u), 
x ■ p Jo 

(7r(p)|M(:r)27 5 d(0)|0> = Jl^L- /' ^^Jm), 

m u + rrid Jo 

(7T(p)\u(x)a^ 5 d(0)\0) = i{p,x u -p u x,) . Um " f'due^M, (3.9) 

6{m u + m d ) Jo 

where ip n is the leading twist 2 distribution amplitude. ip p and ip a are twist 3 distribution 
amplitudes, while g\ and g 2 are of twist 4. Furthermore, we would also include three-particle 
distribution functions [10-12], 
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(n{p)\u(x)a^ 5 F aP (ux)d(0)\0) = -f 3n [{p - PpP^av) - (PaPvgpv 

(7r(p)\u(x)-f^ 5 F a p(ujx)d{0)\0) = if n \pp(g ait - - p a (g^^ - 



x / Va tV± (a t )e^ + ^> x + tf n ^( Pa x^-p^x a ) 
J p ■ x 



■^aPfx * 



Pa{gp„ - — -)] 



(7r(p)\u(x)>y fl F al3 (ujx)d(0)\0} = -f n [pp{g a ^ 

p ■ x p ■ X 

H^+uwfrx _ f Pa_( PaXp - PpXa ) 

p ■ X 

x / Va^ia^^^, 



(3.10) 



where Voti = daida 2 da 3 5(l — a 1 — a 2 — ot 3 ) and F a/3 = \e a p pc F pc . ip 3n is a distribution 
amplitude of twist 3, while <pj_, tp», tp± and <p» are of twist 4. 

After contracting the effective heavy quark fields Qti x i)Qti x 2) m t° propagator 
P + / °° dt5(xi — x 2 — vt), the correlation function (3.8) becomes 



i J d 4 x J d A y dl dte i ^- m ^- x {n{p)\u(x)YP+5{x - y - vl) 

df y) - v a v p df y) d p [y) - &( y) A a {y) - A a (y)d? y) + v a v df y) A^y) 
+v a v p A«(y)dP y) + l -a aP F^{y)]P + 5{y - vth 5 d(0)\0) 



i 
x 



(3.11) 



with &fo = d/(dy a ), which includes only terms related to the two-particle and three-particle 
distribution functions (3.9) and (3.10), while other terms are not considered in this paper. 

Then one can represent the matrix elements by the distribution functions defined in 
Eqs.(3.9) and (3.10). In doing this we choose to work in the fixed-point gauge x p A p (x) = 0, 
where one has 

A p (x) = x v \ dwuFv^iux). 



Using Eqs. (3.5), (3.9) and (3.10), we finally get 



jt^/tt [36/i^(w - l)f f ip p + 18/i^w - 1) if' p - 6^171 J(p a 

+^{u - l)£(10f 2 - 2ml) V ' a + 3^(« - 1) 2 ^V: + 36m^ 2 
-(u - 1)(60£ 2 - 12m 2 )^ - 18(« - lfl 2 g' 2 ' 



+ 



2? 



tzU / da 3 — 
" Jo «3 



2£ (<p ± - <p ± ) + m 7r (^ ± + ^||) 



ci2=l— a\— Qs 



- ai= \— az — ^ 
ai=l-f ' 



(3.12) 



12? 



u 36( M - + 18(« - irrex + 2^e(r + 3m> CT 



-^{u - - 2mi)^ - 3^(u - l) 2 ^Va + 36(C 2 + 3m 2 )^ 
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-108(« - " 18(« - l) 2 fV" - me + 5ml)g 2 

+(u - 1)(72# 2 - 12ml)g' 2 + 18(« - 1) 2 M , 

—kf* [ da 3 I du f dp—{2£ 2 + ml)\ip ± + ip^ 

£ d JO JO JO 0:3 L ll Ja 1 =l-Mpa 3 -^.,Q!2=(Mp-l)o3 + 25 

/ rf «3— [[2y± - I"' ) 7 2 S (3.13) 

where £ 2 = £ 2 — m 2 , and for arbitrary function f(x,y) we use the symbols [f(x,y)] x=Xo = 
/(xo) and [/(#, y)]^=x2 = f{ x i->v) ~ f( x 2,y)- For the detailed techniques used in evaluating 
Eqs.(3.12)-(3.13), we refer to Ref. [4]. 



IV. NUMERICAL ANALYSIS 

With the scale dependence written explicitly, the distribution amplitudes read [2,10-13] 

q -1 r 

<p v (u, pt) = 6u(l -u){l + -a 2 (pt)[5(2u - l) 2 - 1] + — a 4 (/x)[21(2« - l) 4 
-14(2 M -1) 2 + 1]}, 

tp p {u,n) = 1 + ]-B 2 {n)[?,{2u - If - 1] + i J B 4 ( / u)[35(2 M - l) 4 - 30(2« - l) 2 + 3], 
2 o 

tp^n) = 6«(1 - «){l + |c 2 (Ai)[5(2« - l) 2 - 1] + yC 4 (Ai)[21(2« - l) 4 
-14(2«-1) 2 + 1]}, 

^(«, n) = ^ 2 (/V (1 - uf + ^e(fi)5 2 (fi)[u(l -u)(2+ 13«(1 - «)) 

6 

+10w 3 logw(2 - 3m + -m 2 ) + 10(1 - w) 3 log((l - u){2 - 3(1 - u) 

5 

+^(l-«) 2 ))], 
^K/i) = y<5 2 (/i)w(l -u)(2u- 1), 

V? 37r (ai, //) = 360aia 2 a 2 [l + ^i,o(At)(7a 3 - 3) + u 2fi (fi)(2 - 4aia 2 - 8a 3 + 8a 2 ) 

+a;i j i(//)(3aia 2 - 2a 3 + 3a 2 )], 
</?±(«i, /i) = 30<5 2 ( / u)(ai - a 2 )a 2 [i + 2e(^)(l - 2a 3 )], 
(p\\(ai,n) = 1205 2 (fi)e(fi)(ai - a 2 )aia 2 a 3 , 
<P±.{<Xi,li) = 305 2 (/i)a 2 (l - a 3 )[i + 2e(/i)(l - 2a 3 )], 

Cp\\(ai,n) = -120<5 2 ( / u)aia 2 a 3 [i +e(^)(l - 3a 3 )]. (4.1) 

The asymptotic form of these functions and the renormalization scale dependence are given 
by perturbative QCD [11,14]. In the current study, \i should be a typical scale of B — > ixlv 
decay. We set \i to the typical virtuality of the bottom quark, 
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fi b = y/m 2 B - m 2 b « 2.4GeV, (4.2) 

which is the same as that used in Refs. [2,10]. Accordingly, the parameters appearing in the 
distribution amplitudes (4.1) are taken as [2,10] 

a 2 (ii b ) = 0.35, a 4 (/V) = 0.18, B 2 (fi b ) = 0.29, B A (fi b ) = 0.58, 
C 2 (fJib) = 0-059, C 4 (fx b ) = 0.034, 5 2 (fx b ) = 0.17GeV 2 , e(fi b ) = 0.36. (4.3) 

Other meson quantities needed in the sum rule study have been studied via sum rules and 
other approaches. Here we use ^ = 2.02GeV, / w = 0.132GeV [2,10], A B = 0.53GeV and 
F = 0.30GeV 3/2 [8]. 

With these quantities, the variation of 1/rriQ order wave functions with respect to the 
Borel parameter T is shown in Fig.l, where v ■ p = 2.64 GeV is fixed. Appropriate values 
of the threshold s should be determined by the stability of the wave functions. Due to 
the general criterion that both the higher nonperturbative corrections and the contributions 
from excited and continuum states should not exceed 30%, we focus on the region around 
T = 1.5GeV. In this region, as shown in Fig.l, there is no corresponding value of threshold 
which makes the curves of wave funtions very flat. We choose so ~ l-8GeV and s ~ 0.8GeV 
respectively, where the wave functions A 1 and B\ become relatively stable, though the 
stability here is not as optimistic as that for the leading order wave functions A and B [4] . 

Using the thresholds given above, the wave functions are obtained straightforwardly 
from Eqs.(3.6), (3.12) and (3.13). And the form factors f+{q 2 ) and fo{q 2 ) can be easily 
evaluated from the relations (2.2) and (2.19). However, it is known that the light cone 
expansion and the sum rule method would break down as q 2 approaches near ml [2]. At 
large q 2 (or very small y), the results directly evaluated from the sum rule (3.6) turn out 
to be not trustworthy. In order to probe the dynamics of the whole kinematically allowed 
region, we have to extrapolate the sum rule results in the small q 2 region to large q 2 region 
in an appropriate way. For this purpose, for f+{q 2 ) we use the sum rule results at small 
momentum transfer, whereas at large q 2 we assume the single pole approximation [2] 

where rag. = 5.325GeV, /#* ~ 0.16GeV and (]b*Btt ~ 29 [2] would be used. Furthermore, 
it is convenient to parametrize the form factors as 

F(g2) = l-a F q 2 /m 2 B \b F q^/m%' (45) 

where F can be / + and f Q , respectively. We fit the parameters a/ + and bf + by using the sum 
rules and Eq.(4.4), i.e., requiring Eq.(4.5) approach the sum rule results at q 2 < 15GeV 2 but 
be compatible with Eq.(4.4) at q 2 > 15GeV 2 . / is not important in extracting \ V ub \, because 
when the lepton masses are neglected, the B — > ttIu decay width depends only on /+. So we 
use the low q 2 (q 2 < 15GeV 2 ) region results from LCSR to fix the relevant parameters for 
/o in Eq.(4.5). Since /o is nearly constant in q 2 [15], one may expect this parameterization 
yield reasonable approximation for /o in the whole kinematical region. 

In this way we get the fitted parameters in table 1. Then the variation of /+ and f with 
respect to the momentum transfer q 2 are plotted in Fig.2. In table 2 we present values of 
wave functions at some kinematical points of momentum transfer. 
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F(0) 


ap 


7 

b F 




/+ 


0.35 


1.31 


0.35 


LO 




0.38 


1.19 


0.25 


NLO 


fo 


0.35 


0.60 


-0.19 


LO 




0.38 


0.71 


0.41 


NLO 



Table 1. Results of LCSR calculations up to leading (LO) and next leading 
order (NLO) in HQEFT. The leading order results are obtained in Ref. [4]. 
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0.51 


0.72 


1.16 


2.50 


HQEFT-NLO-LCSR 


0.35 


0.48 


0.71 


1.18 


2.70 


HQEFT-LO-LCSR [4] 


0.42 


0.54 


0.77 






HQET-NLO-LCSR [16] 


0.36 


0.47 


0.68 






HQET-LO-LCSR [16] 


0.28 










QCD-LCSR [3] 


0.27 










QCD-LCSR [2 


0.24 










QCD-SR [15] 


fo 


0.38 


0.44 


0.50 


0.54 


0.55 


HQEFT-NLO-LCSR 


0.35 


0.41 


0.50 


0.66 


1.02 


HQEFT-LO-LCSR [17] 



Table 2. Values of the form factors /+ and fo at some kinematical points. 
Neglecting lepton mass, the differential decay width of B — > nlis is given by 

^ 3/2 [/+(? 2 )] 2 



dT Gp\V u i\ 2 , 2 



dq 



24tt 3 



(4.6) 



with E w = (m 2 B + m 2 — q 2 )/(2m B ). Now using the branching ratio Br(B° — > n l + ui) 
(1.8 ± 0.6) x 10~ 4 and the lifetime r B o = 1.542 ± 0.016 ps [18], one gets from Eq.(4.6) 



\V ub \ = (3.2 ±0.5 ±0.2) x 10~ 3 , 



(4.7) 



where the first (second) error corresponds to the experimental (theoretical) uncertainty. 
Here the theoretical uncertainty is mainly considered from the threshold effects. Besides 
this, theoretical uncertainty may arise also from using the single pole approximation to fit 
the form factors at small recoil. To take this into account effectively, if we let the couplings 
in Eq.(4.4) change in the regions /#* = 0.16 ± 0.03GeV and Qb*b-k = 29 ± 3 [2], we get a 
more conservative value for \V u b\\ 

\V ub \ = (3.2 ±0.5 ±0.4) x 10~ 3 . (4.8) 
Eq.(4.7) and (4.8) can be compared with the average obtained by CLEO [18], 

\V ub \ = (3.2518;! ± 0.55) x 10" 3 , (4.9) 
where the uncertainties are statistical and experimental ones. 
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It is also interesting to compare our results with the relations predicted by the large 
energy effective theory (LEET) [19]: 




fVo(. 2 ). 

7T H 



(4.10) 



m\ + ml 



In Fig. 3 we see that the form factor ratio f+/fo derived from our calculations in HQEFT is 
compatible with that predicted by the LEET relation (4.10). They agree well with each other 
at regions not very far from the maximum recoil point. As a comparison, both numerical 
results of the 1/m order contribution in this manuscript and that in Ref. [16] are not large, 
but the former turns out to be a little smaller than the latter. This may be attributed to 
the fact that we have considered the antiquark contributions and also the input values for 
the involved parameters are slightly different. 



We have studied the 1/toq order corrections to B — > ixlv decay in HQEFT, while most 
formulae and discussions here can also be generally applied to other heavy-to-light pseu- 
doscalar meson decays. In HQEFT, 1/mg order corrections from the effective current and 
from effective Lagrangian are given by the same operator forms. Consequently the indepen- 
dent wave functions are reduced in HQEFT in comparison with the ones in HQET. These 
1/mg order contributions have been calculated using light cone sum rule with considering 
two-particle and three-particle distribution amplitudes up to twist 4. Numerically, the 1/mQ 
order wave functions give only about 10% correction to the transition form factor / + . This 
correction indicates a slightly smaller value of the CKM matrix element \V u b\. As to the 
form factor ratio /+/ f , good agreement is observed between the sum rules in HQEFT and 
the LEET at large recoil region. 
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Fig.l. Variation of 1/rriQ order wave functions A\ ((a)) and B\ ((b)) 
with respect to the Borel parameter T at the momentum transfer q 2 = 
OGeV 2 (or £ = v ■ p = 2.64GeV). The dashed, solid and dotted curves in 
(a) correspond to the thresholds so =1-5, 1.8 and 2.1 GeV respectively, 
while those in (b) are derived at sq =0.3, 0.8 and 1.1 GeV respectively. 
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Fig. 2. Form factors /+ ((a)) and fo ((b)) obtained from light cone 
sum rules in HQEFT. The dashed curves are the leading order results 
in HQEFT [4,17], while the solid curves are the results with including 
1/rriQ order correction. 
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Fig. 3. Comparison of the ratio f+(q 2 )/fo(q 2 )- The dashed curve is 
obtained from the LEET relation (4.10), and the dotted curve is from 
HQEFT sum rule with including only the leading order wave functions 
[4,17], while the solid curve is the result of this paper, which includs also 
the 1/rriQ order contributions. 
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